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•  Generalized  3D  computational  model  of  an  electric  double  layer  supercapacitor. 

•  3D  microstructural  aspects  do  not  have  a  significant  impact  on  the  performance. 

•  Specific  capacitance,  ionic  conductivity,  and  tortuosity  are  critical. 

•  State-of-the-art  numerical  methods  provide  accurate  and  robust  solutions. 
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A  generalized  three  dimensional  computational  model  based  on  unified  formulation  of  electrode 
—electrolyte  system  of  an  electric  double  layer  supercapacitor  has  been  developed.  This  model  accounts 
for  charge  transport  across  the  electrode-electrolyte  system.  It  is  based  on  volume  averaging,  a  widely 
used  technique  in  multiphase  flow  modeling  ([1,2])  and  is  analogous  to  porous  media  theory  employed 
for  electrochemical  systems  [3-5],  A  single-domain  approach  is  considered  in  the  formulation  where 
there  is  no  need  to  model  the  interfacial  boundary  conditions  explicitly  as  done  in  prior  literature  ([6]). 
Spatio-temporal  variations,  anisotropic  physical  properties,  and  upscaled  parameters  from  lower  length- 
scale  simulations  and  experiments  can  be  easily  introduced  in  the  formulation.  Model  complexities  like 
irregular  geometric  configuration,  porous  electrodes,  charge  transport  and  related  performance  char¬ 
acteristics  of  the  supercapacitor  can  be  effectively  captured  in  higher  dimensions.  This  generalized  model 
also  provides  insight  into  the  applicability  of  ID  models  ([6])  and  where  multidimensional  effects  need 
to  be  considered.  A  sensitivity  analysis  is  presented  to  ascertain  the  dependence  of  the  charge  and 
discharge  processes  on  key  model  parameters.  Finally,  application  of  the  formulation  to  non-planar 
supercapacitors  is  presented. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Ultracapacitors  or  supercapacitors  are  charge  storage  devices 
that  operate  on  the  principle  of  electrochemical  double  layer 
(ECDL)  capacitance  wherein,  electrical  energy  can  be  stored  and 
released  by  nanoscale  charge  separation  at  the  interface  between 
the  electrode  and  the  electrolyte.  In  paper  [7]  advances  in  electrode 
materials  for  the  supercapacitor  are  summarized.  There  has  been 
considerable  activity  in  the  recent  years  to  exploit  high  surface  area 
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offered  by  Nanotubes  [8-14]  and  more  recently  Graphene  [15], 
Most  of  the  projected  gains  in  energy  densities  (getting  closer  to 
that  of  current  Li-ion  batteries)  are  due  to  increase  in  capacitance 
through  increase  in  surface  area  and  exploiting  the  nano-pore  and 
ion  interactions.  However,  considerable  development  needs  to  be 
done  at  the  system  level  to  ascertain  the  true  performance  in  a 
practical  supercapacitor  [16],  In  this  work,  we  developed  a 
macroscopic  model  that  can  rapidly  use  the  microscale  properties 
to  assess  the  overall  performance  of  the  supercapacitor  device  and 
possibly  aid  in  transition  of  this  rapid  progress  in  supercapacitor 
electrode  materials  into  high  performance  supercapacitors  without 
lot  of  iterations  at  the  device  level  where  one  needs  to  balance  the 
cathode,  anode,  electrolyte,  and  current  collectors  to  maximize 
energy  and  power  density. 
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Advantages  in  using  the  nanoporous  carbon  materials  with  fine 
tuned  pore  sizes  for  the  corresponding  electrolytes  is  well  estab¬ 
lished.  Especially,  the  cathode  and  anode  with  different  pore  sizes 
that  match  the  anion  and  cation  sizes  could  lead  to  improved  power 
delivery.  Aspects  such  as  high  conductivity  and  surface  area  of 
these  carbon  material  characteristics  are  studied  for  various  engi¬ 
neered  carbon  material  with  tailored  pore-size  distribution  for  high 
capacitance  [10,11,17], 

The  amount  of  energy  stored  in  a  supercapacitor  is  inversely 
proportional  to  the  thickness  of  the  double  layer  and  directly  pro¬ 
portional  to  the  specific  surface  area  of  the  electrodes.  By  using 
nanoporous  electrodes,  supercapacitors  deliver  a  sufficiently  high 
power  density  compared  to  conventional  batteries  and  an 
extremely  high  energy  density  compared  to  conventional  dielectric 
capacitors.  However,  the  energy  density  of  supercapacitor  devices 
are  an  order  of  magnitude  lower  than  current  battery  devices. 

They  are  also  amenable  to  a  variety  of  applications  ranging  from 
energy-smoothing,  regenerative  braking  to  power-source  in  hybrid 
vehicles.  Supercapacitors  have  the  advantage  of  long  operating  life, 
a  wide  thermal  operating  range,  low  weight,  low  maintenance  and 
near  infinite  cyclability  (charge— discharge). 

A  Ragone  plot  (see  Fig.  1)  shows  the  power  density  vs.  energy 
density  for  common  energy  storage/conversion  devices. 

Currently,  ultracapacitors  need  coupling  with  batteries  to  sup¬ 
ply  energy  for  a  prolonged  time  period  [  18  ].  One  of  the  motivations 
for  developing  novel  materials  is  to  increase  their  energy  density 
and  better  simulation  models  can  accelerate  the  adoption  of  the 
most  promising  materials  from  a  full  system  perspective.  In  addi¬ 
tion,  several  new  designs  at  the  system  level  can  be  explored  to 
maximize  the  energy  density  while  retaining  very  high  power 
density  and  life.  This  requires  a  careful  understanding  of  the  physics 
at  the  nanoscale  as  well  as  the  device  scale.  While  there  is  a  sig¬ 
nificant  increase  in  research  at  the  nanoscale  [11,12,19,20],  there  is 
much  work  to  be  done  at  the  device  scale.  Even  recent  publications 
e.g.  Ref.  [6]  still  consider  the  electrode  properties  to  be  isotropic  and 
the  simulations  are  done  in  one-dimension.  With  the  advent  of  new 
graphene  based  electrode  materials  that  are  utilized  for  the  ultra¬ 
capacitor  [21  ],  due  to  its  high  surface  area  such  as  2630  m2  g_1  [22], 
the  accuracy  of  these  assumptions  with  reference  to  truly  nano- 
structured  materials  like  chemically  modified  graphene  and  carbon 
nanotube  solids  have  not  been  addressed.  Also,  most  of  these 
models  assume  the  capacitance  in  the  electrode  pores  (due  to 
ECDL)  can  be  modeled  using  parallel-plate  assumption.  It  has  been 
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Fig.  1.  A  Ragone  chart  showing  the  power  vs.  energy  densities  for  a  number  of  energy 
storage  devices  and  compares  that  to  a  gasoline  engine.  It  can  be  noticed  that  ultra¬ 
capacitors  have  almost  1000  times  the  energy  density  of  dielectric  capacitors  and 
around  10  times  the  power  density  of  conventional  batteries.  Units  in  this  figure  have 
been  modified  to  be  consistent  with  the  journal  requirement. 


shown  recently  that  capacitance  exhibits  a  strong  dependence  with 
pore  size  and  curvature  [11,20],  Lanzi  and  Landou  [23]  show  mac- 
roporous  structure  is  important  in  determining  the  ohmic  and 
mass-transfer  resistance  whereas  in  micropores  these  effects  are 
negligible  because  transport  is  too  localized.  They  also  find  that 
distribution  in  pore  volume  between  macropores  and  micropores 
can  have  a  major  effect  on  porous  electrode  performance.  Increase 
in  porous  material  surface  area  for  different  activated  carbon  ma¬ 
terial  and  its  influence  to  the  capacitor  performance  i.e.,  specific 
capacitance,  one  of  the  important  characteristics  of  the  carbon 
electrode  material,  has  been  established  [24,25], 

In  this  paper,  we  test  the  validity  of  the  first  assumption,  namely, 
the  effect  of  anisotropy  in  electrode  properties  on  the  capacitance, 
charge  and  discharge  characteristic  of  a  ultracapacitor. 

For  our  model,  we  consider  an  ultracapacitor  unit-cell  as  shown 
in  Fig.  2  comprising  of  two  nanoporous  electrodes  isolated  from 
electrical  contact  through  a  porous  separator.  The  pores  are  filled 
with  an  electrolyte  that  allows  ionic  current  to  flow  through  the 
electrodes  while  preventing  self-discharge.  Current  collectors 
typically  made  of  metal  foils  drain  the  electrical  current  from  each 
electrode.  In  a  full-scale  system,  multitude  of  these  unit-cells  are 
used.  For  our  modeling  purpose,  it  is  sufficient  to  characterize  one 
of  these  unit-cells. 

A  historical  perspective  on  modeling  electrode-electrolyte 
interface  as  circuit  has  been  presented  by  Geddes  [26],  The 
macroscopic  models  are  primarily  based  on  electrical  circuit 
models  or  simple  lD/pseudo-2D  models  [27,28],  V.  SrinivasanandJ. 
W.  Weidner  [29]  developed  analytical  models  for  the  capacitors 
under  constant  current  while  neglecting  the  electrode  kinetics. 
Pillay  and  Newman  [30]  have  studied  the  influence  of  side  reactions 
using  Tafel  kinetics  on  performance  of  EDCL.  A  mathematical  model 
is  proposed  by  Ong  and  Newman  [31  ]  for  double  layer  capacitor  and 
the  predictions  show  slower  charge  build-up  due  to  double  layer 


Fig.  2.  Schematic  of  an  electrochemical  double-layer  supercapacitor.  The  electrodes 
are  typically  made  of  a  nanoporous  material.  The  porous  separator  and  the  electrode 
pores  are  saturated  with  the  electrolyte. 
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capacitance  during  galvanostatic  charge.  Also  their  analysis  indi¬ 
cated  no  difference  to  the  chemical  kinetics  model  on  the  potential 
profiles.  A  new  RC  circuit  model  was  proposed  by  Dunn  and 
Newman  [32]  for  the  double-layer  capacitors  and  maximum 
theoretical  specific  energy  over  a  range  of  discharge  times  for 
constant-current  and  constant-power  discharges  were  predicted 
by  optimizing  the  electrode  thickness,  electrode  porosity,  and  final 
voltage.  Further  modeling  was  done  by  Verbrugge  and  Liu  [6]  to 
include  the  microstructure  and  pore  size  distribution  in  calculating 
the  effective  properties  of  the  porous  electrodes. 

The  modeling  technique  adopted  in  these  works  allows  for 
splitting  the  operator  in  each  domain.  Then  solution  to  these  in¬ 
dividual  operators  are  obtained  by  solving  them  in  segregated 
fashion.  Though  solutions  are  obtained  quickly,  these  segregated 
solution  techniques  give  rise  to  numerical  instabilities  and 
convergence  issues.  The  main  problem  might  be  the  fact  that  the 
strongly  elliptic  equations  have  been  solved  as  parabolic/hyperbolic 
equations  by  segmenting  the  domain  and  specifying  hand-shake 
boundaries  at  the  interfaces.  This  leads  to  ill-posed  numerical 
scheme  where  iterative  techniques  may  not  converge  and  solutions 
are  not  assured  in  all  possible  problem  setups. 

These  models  are  limited  in  their  ability  to  capture  spatial  var¬ 
iations  in  permeability  or  conductivity  or  able  to  handle  multidi¬ 
mensional  structure  of  the  electrode  material.  Any  variations  in  the 
porosity  of  the  solid  material  or  anisotropy  of  the  physical  prop¬ 
erties  are  not  accounted.  All  previous  models  have  assumed  a  single 
value  of  electrode  porosity  and  have  not  considered  the  effect  of 
varying  pore  structure,  electrode  composition  on  electrochemical 
capacitor  performance  [33],  Also,  the  temporal  changes  of  the 
transport  properties  are  neglected  in  the  simulations. 

1.1.  Scope  of  present  work 

In  this  paper  we  propose  a  unified  multi-dimensional  formu¬ 
lation  across  the  electrode-electrolyte-electrode  system  using  a 
volume  averaging  approach  typical  of  multiphase  flow  formulation 
([1,2])  and  is  similar  to  the  porous  media  theory  often  used  for 
modeling  batteries  [3—5],  Following  are  some  of  its  advantages 
over  the  ID  models: 

•  General  multi-dimensional  formulation  allowing  for  realistic  2D 
and  3D  effects 

•  Anisotropic  electric  and  heat  conductance,  capacitance,  chem¬ 
istry  etc.  can  be  specified  as  a  function  of  space  and  time 

•  This  formulation  is  also  suitable  for  any  upscaling  from  lower 
length-  and  time-scale  methods  as  the  properties  can  be 
determined  locally  (at  the  discretization  level) 

•  Because  of  the  single-domain  approach,  complex  geometries  are 
naturally  incorporated.  The  same  formulation  applies  to  ID,  2D 
and  3D  geometries  irrespective  of  the  complexities  associated 
with  dimensionality  as  well  as  with  electrode/electrolyte  spatial 
arrangements 

•  Unlike  previously  used  methods  for  supercapacitors  based  on  a 
segregated  formulation  with  intermediate  boundaries,  by  hav¬ 
ing  a  unified,  single-domain  approach,  stability  and  conver¬ 
gence  of  the  numerical  algorithm  is  fundamentally  guaranteed 

•  Implicit  coupled  solutions  can  be  obtained  with  improved  nu¬ 
merical  efficiency  and  can  be  scaled  to  run  of  massively  parallel 
computers  as  needed 

We  first  validate  this  unified  approach  by  comparing  with  ex¬ 
periments  reported  in  the  literature  to  verify  the  accuracy  and  ef¬ 
ficiency  of  the  method.  Furthermore,  we  have  tested  some  of  the 
assumptions  made  currently  in  terms  of  isotropy,  homogeneity  etc., 
to  understand  the  sensitivity  to  the  overall  performance  of  the 


supercapacitor  but  also  to  the  internal  fields  such  as  local  potential 
and  current  within  the  cell.  We  also  applied  this  method  to  a  2D 
supercapacitor  configuration  to  study  the  effects  of  non-planar 
arrangements. 

To  arrive  at  the  governing  equations  at  the  device  scale,  we  use 
the  volume  averaging  approach  wherein,  the  electrode/separator 
and  the  electrolyte  are  considered  as  inter-penetrating  continua. 
Individual  conservation  equations  for  each  phase  are  averaged  over 
a  volume  element  that  is  larger  than  the  pore  scale  and  smaller  than 
the  device  scale  (here,  the  electrode  thickness).  Clearly,  this  re¬ 
quires  that  the  microscopic  features  of  the  electrode  be  more  or  less 
uniform  and  at  a  length  scale  much  smaller  than  the  rest  of  the 
domain.  This  requirement  can  be  violated  for  truly  nanostructured 
materials  e.g.,  in  CNT  solids,  the  pores  can  extend  right  through  the 
electrode.  Thus  averaging  along  the  thickness  of  the  electrode  is  not 
valid.  This  is  an  extreme  case  and  we  shall  assume  that  the  effect  of 
nanoscale  structures  can  be  captured  by  considering  anisotropy  in 
the  electrode  material  properties. 

We  also  assume  the  dilute-solution  theory  for  the  liquid  phase 
in  the  pores  of  the  electrode.  This  assumption  is  validated  since  the 
salt  concentration  in  the  electrolyte  is  small  enough  (<1  M)  and 
that  the  electrolyte  conductivity  is  a  linear  function  of  the  salt 
concentration  in  that  regime  [6],  Furthermore,  little  information  is 
available  on  the  concentration  dependence  of  the  physical  and 
chemical  parameters  of  the  system.  As  more  information  about  the 
parameters  is  available  for  e.g.,  through  ab-initio  simulations,  this 
assumption  has  to  be  revisited. 

This  paper  is  structured  as  follows:  In  the  next  section,  a 
mathematical  model  for  the  ECDL  is  provided  with  a  brief  discus¬ 
sion  on  its  validity.  Following  this,  the  analysis  for  the  case  of  a 
uniform  property,  isotropic,  binary  electrolyte  is  provided.  Results 
are  compared  with  literature.  We  then,  consider  the  case  of 
anisotropy  and  spatial  variability  in  the  matrix  properties  (which 
also  effects  liquid  phase  properties  through  porosity  variations). 
Finally,  we  present  results  from  a  2D  supercapacitor  configuration 
followed  by  conclusions. 

2.  Mathematical  model 

To  derive  the  volume-averaged  mathematical  model,  we  start 
with  the  conservation  of  charge  for  a  infinitesimal  volume  element 
inside  the  supercapacitor: 

|q  +  V-J  =  0  (1) 

where,  q  is  the  charge  contained  in  the  volume  element  and  J  is  the 
current  density  that  acts  as  a  flux  term  on  the  boundaries  of  the 
volume  element.  By  Nernst— Planck  relation,  for  a  given  electrolyte, 
there  are  three  possible  modes  that  can  affect  the  current  density: 
a)  voltage  gradients,  b)  concentration  of  different  species,  and  c) 
external  velocities  (for  species  convection).  In  this  paper,  for 
simplicity,  we  only  consider  supercapacitors  with  no  chemical  re¬ 
actions  or  external  forcing  terms  and  the  framework  can  be  easily 
extendable  to  include  chemical  reactions  and  external  forcing. 
Therefore,  the  current  density  can  be  modeled  according  to  Ohm’s 
law.  Before  we  can  derive  the  supercapacitor  model,  we  need  a  few 
preliminaries: 

Volume  fraction:  The  volume  fraction  of  a  component  is  calcu¬ 
lated  as 


where,  dV  is  the  volume  of  the  infinitesimal  volume  element  and 
d\4  is  the  volume  occupied  by  the  k-th  phase/species  in  that 
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element.  Let  us  now  consider  the  volume  average  of  the  charge 
conservation  equation  Eq.  ( 1 ) 

a(q)  +  (V-J>  =  0, 

where  averaging  operator(g)  =  J  gdQ 
dv 

In  the  L.H.S,  it  is  tacitly  assumed  that  the  volume  element  does 
not  vary  with  time,  i.e.,  there  is  no  dilatation.  Let  us  now  consider 
physical  models  for  the  various  terms  in  the  above  equation. 

Model  for  q:  In  a  supercapacitor,  the  capacitance  arises  from 
charge  separation  across  the  double  layer  that  occurs  at  the  inter¬ 
face  of  the  electrode  and  electrolyte.  Let  the  local  capacitance  in  the 
volume  element  be  C  and  the  local  voltage  across  the  double  layer 
be  V.  In  this  volume  averaged  approach,  we  are  capturing  the  effect 
of  the  double  layer  while  explicitly  not  resolving  the  interface  be¬ 
tween  the  electrode  and  electrolyte.  We  have  local  charge  density 
as  q  =  CV.  On  taking  the  volume  average,  we  get 
(q)  =  (C)(V)  +  (C'V1).  Where,  C'  and  v'  are  the  fluctuations  in  the 
local  capacitance  and  voltage  due  to  spatial  variations  of  the  elec¬ 
trode  microstructure  and  properties  below  the  resolution  of  the 
volume  element.  Currently,  only  the  averaged  capacitance  of  a 
supercapacitor  is  measured  and  these  spatial  variations  are  not 
available.  One  could  perform  simulations  analogous  to  DNS  (Direct 
Numerical  Simulations)  for  batteries  [34],  where  the  supercapacitor 
morphology  along  with  the  double  layer  can  be  explicitly  resolved 
to  provide  these  closures.  For  now,  we  assume  that  C  =  0  and  this 
additional  term  can  be  included  once  we  have  more  data.  This 
yields: 

(i q )  =  CV,  where  C  =  <C),V  =  (v^  (4) 

Model  for  J:  The  Ohm’s  law  for  a  given  component  (electrode  or 
electrolyte)  can  be  written  to  be  of  the  following  form 


term.  For  a  full-cell  simulation,  one  needs  to  formulate  in  terms  of 
the  potential  in  the  neutral  phase  of  the  electrolyte  and  electrode 
materials.  We  note  that  conservation  of  current  is  ensured  i.e., 
i'i  +  fa  =  0  and  no  species  gradient  is  assumed  which  is  true  for  the 
region  outside  of  electric  double  layer.  These  governing  equations 
can  be  written  as  follows: 


Clt{Vl~Vs)  *  V<k,e,Wi)  (8) 

C^M-V,)  =  V-(ksesVVs)  (9) 

This  set  of  equations  model  any  change  of  potential  in  electrode 
phase  and  a  corresponding  automatic  response  of  the  potential  in 
electrolyte  phase  to  ensure  charge  neutrality.  The  multi-domain 
specific  equations  used  in  Ref.  [6]  reduce  to  this  form  and  both 
formulations  are  consistent.  The  biggest  advantage  of  the  current 
formulation  is  that  these  set  of  equations  are  applicable  over  the 
entire  supercapacitor  and  the  local  volume  fraction  and  corre¬ 
sponding  phasic  properties  ensures  the  consistency  between  these 
two  different  formulations.  We  can  therefore  only  specify  the 
boundary  conditions  at  the  current  collectors  and  treat  complex  2D 
and  3D  configurations  easily  as  we  do  not  have  to  explicitly  handle 
internal  domain  boundaries. 


2.1.  Numerical  solution  procedure 

In  this  section  we  briefly  describe  the  solution  procedure 
employed  in  solving  the  system  of  Eqs.  (8)  and  (9)  for  charge 
conservation  in  electrochemical  double  layer  super-capacitor. 

These  equations  are  discretized  in  time  using  first-order  BDF 
(backward  differentiation  formula)  and  in  space  using  standard 
finite  element  method.  For  the  discretization  based  on  finite 
element  formulation  we  begin  with  the  system  of  equations, 


Jot  =  -kaVVa  (5) 

where,  Ja  is  the  current  density,  Va  is  the  local  potential,  and  ka  is 
electrical  conductivity  within  the  material  component  a.  By  volume 
averaging  the  above  equation  and  utilizing  the  definition  of  volume 
fraction,  we  get: 


V-(k,e,  W,)j  ,n  Q 
C&(Vs-Vi)  -  V-(ksesWs)J 


with  boundary  conditions, 
for  galvanostatic  charge: 


(10) 


(J)  =  -(kafa)  VVa  +  (V  W£)  (6) 

where  the  fluctuation  term  again  represents  the  unresolved  fluc¬ 
tuations  due  to  the  variations  in  microstructure  and  properties 
below  the  volume  element  that  we  do  not  explicitly  model.  Though 
this  term  could  be  of  considerable  importance  to  understand  local 
variations,  due  to  lack  of  knowledge  about  the  microstructure  of  the 
electrodes  or  how  the  properties  vary  at  those  scales,  we  neglect 
this  term.  Since,  we  are  interested  in  the  mean  response  of  the 
supercapacitor  at  the  coarse  scales  and  the  dominant  (diffusion) 
operator  is  linear,  this  is  not  a  bad  assumption  in  terms  of  capturing 
overall  response  of  the  system.  Here  we  wanted  to  have  a  general 
derivation  so  that  we  can  include  additional  details  when  available 
and  also  provide  a  way  to  link  to  lower  length-scale  simulations  or 
experiments  that  resolve  the  small  scale  variations. 

This  brings  us  to  the  model  equation 

C^V  =  f-((kaea)Wa)  (7) 


°seJTt „  on  A 

=  io  on  r2 

for  potentiostatic  discharge: 

Vs  =  —  Vq/2  on  r ! 

Vs  =  Vq/2  on  r2 


(11) 


(12) 


where  r1TjT2  =  T  =  dfi  and  Tj  nr2  =0 

We  construct  an  integral  form  based  on  Galerkin  Method  with 
weak  form  by  transferring  one  order  of  differentiation  from 
V  =  {Vi,Vs}  to  test  function  w, 

/  (C^8T^)wdQ+  /(fcie/wiVw)dQ  =  w/</erW,  |r 
°  (13) 

J  (C~^T^)wdfi+  J  (kscsVVsVw)dQ  =  wksesVVs  Ir 

n  q 


The  above  equation  is  valid  over  each  component  and  assumes  C 
does  not  vary  with  time.  If  C  varies  with  time,  the  above  formula¬ 
tion  will  be  slightly  modified  with  C  moving  into  the  time  derivative 


This  constitutes  the  weak  form  in  Qx.  Let  Qx  —  U  eQx  be  the 
discretization  of  Qx  in  which  Qx  is  an  element  e.  Now  consider  the 
element  approximation  for  {V}  such  that 
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veh(x,t)  =  jrNi(x)V?(t )  =  [N(x)]\7e(t) 

Using  this  Eq.  (14)  we  can  write 


(14) 


^  =  E  8ia?Vc(t)  =  aJV 

'n1  (15) 

ENiWf^v' 

i  —  1 

Substituting  these  relations  into  Eq.  (13)  we  get, 

/  o4K*r-*o<m.+  /*"$  (^)dn« 

=  tyXe+i^qVV^  -  Nj(xe)k,ei (16) 

J  CNM^.-Vy K>.+ 

=  Nj(xe+1)fesesWs|x  ^  -  Ni(Xe)/csesWsj3f  (17) 

This  can  be  written  as, 

c[/fc]  [Me]{l7eJ  +  [De]{l/e>  =  {Pe}  (18) 


where  Ik  is  the  block  matrix  for  each  component  {Vf,  Vse}  given  by, 


For  discretization  in  time,  we  use  standard  first  order  BDF 
scheme  on  {V  }  which  results  into  following  algebraic  system, 


(20) 

The  strongly  coupled  formulation  is  implemented  within 
Advanced  Multi-Physics  (AMP)  software  framework  [35],  We  use 
the  standard  interface  to  the  Trilinos  direct  sparse  linear  solver  [36] 
through  AMP  to  solve  the  coupled  algebraic  system  (20)  for 
{W}  =  {Vf,Vt}. 

3.  Results 

In  this  section,  we  use  a  commercial  pseudocapacitor  (see 
Table  1  for  baseline  properties)  as  a  demonstration  vehicle  to  pre¬ 
sent  a  number  of  studies  that  employ  the  formulation  detailed 
above  and  highlight  the  flexibility  and  applicability  of  the  proposed 
model.  First,  we  validate  the  charging  and  discharging  profiles  by 


Table  1 

Properties  and  parameters  for  the  simulations. 


Parameter  Value 

Electrode  width,  We  50 

Separator  width,  W5  25 

Electrode  height,  I  25 

Specific  capacitance,  aC  86 

Void  volume  fraction  in  electrodes,  es  0.67 

Void  volume  fraction  in  separator,  c(  0.6 

Electrolyte  conductivity,  k  0.067 

Solid  phase  conductivity,  a  52.1 

Tortuosity  for  liquid  phase  in  carbon  electrode,  r,  2.3 

Tortuosity  for  liquid  phase  in  separator,  n  1.29 


Units 


comparing  to  the  experimental  results.  We  will  also  present  the 
grid  refinement  and  convergence  studies  to  demonstrate  the  ac¬ 
curacy  of  the  solutions.  Subsequent  to  these  studies  we  will  present 
a  set  of  parametric  studies  to  show  the  sensitivity  of  the  results  to 
different  properties,  anisotropy  of  both  physical  and  geometric 
properties,  and  effects  of  non-planar  architectures. 

The  simulations  are  all  carried  out  with  a  current  of  100  A  for  a 
variety  of  charging  and  discharging  voltages.  The  numerical  studies 
are  divided  into  following  sections. 


3.1.  Experimental  validation 

For  model  validation,  we  use  the  data  from  the  3500  F  SAFT 
America  supercapacitor  (Maximum  voltage  2.8  V)  that  was  exten¬ 
sively  studied  earlier  by  Verbrugge  and  Liu  [6],  This  single  cell 
supercapacitor  consists  of  two  carbon  electrodes  and  a  separator 
preventing  physical  contact  between  them.  The  physical  properties 
used  for  the  simulation  are  mostly  derived  from  Ref.  [6]  and  can  be 
found  in  Table  2. 


3.1.1.  Comparison  against  the  results  from  Verbrugge  and  Liu  [6] 

In  Fig.  3,  we  compare  the  simulated  supercapacitor  voltage 
against  the  experimental  results  published  in  Ref.  [6]  for  constant 
current  charging  from  1.6  V  to  2.2  V.  For  the  first  18  s  of  the 
simulation  during  this  charge  phase,  the  cell  voltage  increased  from 
1.6  V  to  2.2  V.  The  charging  time  corresponding  to  the  threshold 
voltage  of  2.2  V  matches  very  well.  At  this  point,  the  boundary 
conditions  are  switched  to  constant  voltage  discharge  at  1.4  V.  In 
Fig.  4,  we  see  a  comparison  of  the  current  from  the  simulations  to 
that  of  the  experiment.  The  rapid  initial  rapid  discharge  followed  by 
slow  discharge  asymptotically  approaching  zero  current  are  accu¬ 
rately  captured  by  the  simulations  and  are  in  very  good  agreement 
with  the  experiments.  This  analysis  gives  us  confidence  that  the 
current  formulation  captures  the  macroscopic  response  of  the 
supercapacitor  very  well  and  we  can  use  this  model  to  further 
understand  the  sensitivities  to  grid  refinement,  model  parameters, 
etc. 


Table  2 

Boundary  conditions. 

Left  current  collector  (x  =  0) 
Galvanostatic  charging 
6V,l6x  =  0 
Vs  =  0 

Potentiostatic  discharging 

SVi/dx  =  0 

Vs  =  -Vdischarge/2 


Right  current  collector  ( x  =  I) 


6V,/6x  =  0 
SVs/Sx  =  iol(ae) 

5V,l5x  =  0 
Vs  =  Vdisch^e/2 
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Fig.  3.  Comparison  of  cell  potential  with  the  Verbrugge  and  Liu  [6]  results.  The  experiments  are  denoted  by  filled  circles  and  simulations  are  +  symbol. 
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Fig.  4.  Comparison  of  discharging  current  with  the  Verbrugge  and  Liu  [6]  results.  The  experiments  are  denoted  by  filled  circles  and  simulations  are  +  symbol. 


3.1.2.  Grid  refinement  and  convergence 

In  Fig.  5,  we  show  a  mesh  refinement  study  conducted  to  eval¬ 
uate  the  convergence  of  the  solution.  For  this  study  we  have  used 
three  different  meshes  with  1200  and  2500  elements  across  the 
domain.  The  residual  in  the  solution  is  of  the  O(10-8)  during  the 
entire  transient  analysis.  This  shows  that  we  have  obtained  grid 
independent  solutions  and  any  further  refinements  would  not 
change  the  solution  significantly. 

The  supercapacitor  problem  we  are  modeling  is  essential  1 D  and 
we  are  modeling  it  as  a  2D  system  with  few  elements  in  the 
transverse  direction.  For  generality  one  could  also  model  it  as  a  3D 
system.  In  Fig.  6,  we  compare  the  cell  potential  and  current  density 
between  3D  and  2D  domains.  As  expected,  the  solutions  are  nearly 
identical  as  there  are  no  gradients  imposed  in  the  third  dimension. 
This  verifies  that  the  formulation  for  higher  dimensions  yields  the 
same  result  as  2D  and  ID  if  the  inherent  problem  is  of  a  lower 
dimension.  This  verification  test  ensures  that  both  the  formulation 
and  the  implementation  are  consistent  to  the  physical  processes  we 
are  modeling  in  the  device. 


3.1.3.  Spatial  variation  of  the  potential  and  current  density 

A  series  of  charge  and  discharge  conditions  are  imposed  to 
evaluate  the  solution  characteristics  of  the  supercapacitor.  Initially, 
a  potentiostatic  charge  is  imposed  to  maintain  a  cell  potential  of 
1.6  V  and  to  equilibrate  any  initial  gradients  in  the  cell.  At  this 
stage,  a  galvanostatic  charge  of  100  A  i.e,  current  density 
64.93  Am-2  is  imposed  for  charging  the  cell.  This  charging  is 
continued  until  the  cell  potential  reaches  2.2  V.  The  cell  potential 
(difference  across  the  cell)  varies  in  time  from  1.6  to  2.2  V  but 
there  are  internal  voltage  gradients  due  to  finite  transport  co¬ 
efficients.  Fig.  7  shows  these  variations  along  the  width  of  super¬ 
capacitor  and  small  gradients  in  voltage  are  established  —  these 
gradients  are  very  minimal  compared  to  overall  cell  voltage.  After 
the  cell  reaches  2.2  V,  the  cell  is  subjected  to  a  potentiostatic 
discharge  at  1.4  V.  As  can  be  seen  in  Fig.  8,  the  gradients  of  the 
potential  within  the  cell  can  be  quite  severe  because  of  high  cur¬ 
rent  discharge  during  the  switch  from  galvanostatic  charge  to 
potentiostatic  discharge.  Usually  gradients  are  detrimental  in 
terms  of  life  and  efficiency  (under  utilizing  electrodes  and  losses) 


Fig.  5.  Mesh  convergence  study. 
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Fig.  6.  Comparison  of  the  charge/discharge  profile  from  3D  to  2D  domains. 


and  simulations  such  as  these  can  identify  operating  conditions  to 
minimize  spatial  gradients. 

3.2.  Sensitivity  analysis  to  physical  properties 

In  this  section,  we  look  at  the  sensitivity  of  the  supercapacitor 
performance  to  variations  in  electronic  conductivity  in  the  elec¬ 
trode  (cr),  ionic  conductivity  in  the  electrolyte  (k),  specific  capaci¬ 
tance  (aC),  and  Bruggeman’s  coefficient  (tortuosity  of  the  pore 
structure)  as  these  are  properties  one  can  dial  in  using  various 
advances  in  materials  science.  The  purpose  of  such  a  parametric 
analysis  is  many  folds:  a)  to  understand  what  properties  are  to  be 
measured  accurately  to  get  accurate  simulation  results,  b)  narrow 


down  fundamental  exploration  to  concentrate  on  variations  that 
will  have  maximum  impact  on  the  overall  system  performance  and 
c)  to  understand  the  relative  interaction  between  the  various 
competing  properties.  Below,  we  fix  the  void  fraction  in  electrodes 
(es),  void  fraction  in  separator  (e/)  corresponding  to  the  base  case  we 
detailed  before  and  vary  one  of  the  properties  (cr,  k,  or  aC)  relative  to 
the  base  case.  One  could  perform  second  order  variations  by 
changing  two  or  more  properties  at  the  same  time  but  we  have  not 
done  that  in  this  study. 

Porous  carbon  with  high  surface  area  are  prepared  by  various 
manufacturing  processes  [37,38],  The  effective  electronic  con¬ 
ductivity  of  these  materials  hugely  depends  on  the  pore  structure 
and  compositions  [39],  For  the  various  configurations  of 


Fig.  7.  Spatial  variation  of  the  local  electrode  potential  during  the  constant  current  charge  : 


A  from  1.6  to  2.2  V. 
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mesoporous  carbon  the  electronic  conductivity  ranges  can  be 
found  in  Ref.  [40],  In  Table  3,  we  list  test  cases  with  varying 
electronic  conductivity  and  study  its  effect  on  the  supercapacitor 
performance.  As  can  be  seen  from  these  simulations  (Fig.  9),  any 
conductivity  greater  than  1.0  S  m-1  gives  almost  identical  results 
and  the  simulations  indicate  that  the  conductivity  has  to  decrease 
to  0.1  S  m-1  before  one  can  see  adverse  affects  and  solid  phase 
conductivity  becomes  rate  limiting.  From  a  practical  perspective 
most  carbons  have  electronic  conductivity  greater  than  1.0  S  m-1 
and  any  further  improvements  in  electronic  conductivity  would 
not  have  any  improvement  in  overall  supercapacitor  performance 
for  typical  applications. 


Table  3 

List  of  cases  for  parametric  study  of  a. 


Parameter  Base  case  Test  case-1 


a  (S  m-1)  52.1  1.0 

K  (S  m  ')  0.067  0.067 

aC  (F  m-3)  86  x  106  86  x  106 

es  0.67  0.67 

e,  0.6  0.6 


100.0 
0.067 
86  x  106 


We  varied  the  ionic  conductivity  of  the  electrolyte  as  part  of  the 
next  study.  Ionic  conductivity  is  usually  the  rate  limiting  step  in 
batteries  (for  e.g.,  solid  state  vs.  liquid  electrolyte)  and  we  wanted 
to  explore  the  dependence  in  supercapacitors.  In  Table  4,  we  list 
test  cases  with  varying  ionic  conductivity  and  study  its  effect  on  the 
supercapacitor  performance.  As  in  the  case  of  batteries,  the  ionic 
conductivity  turns  out  to  be  the  limiting  factor  for  charge  transport, 
any  minimal  variation  in  the  value  shows  a  significant  change  in  the 
solution  and  charge/discharge  rates.  Fig.  10  shows  that  energy 
stored  and  total  capacity  of  supercapacitor  is  drastically  low  at 
lower  values  of  ionic  conductivities.  During  galvanostatic  charge, 
the  cell  reaches  the  cut-off  potential  at  a  fraction  of  time  at  lower 
ionic  conductivity  because  of  the  increased  transport  resistance 
and  during  the  discharge,  the  limited  charge  stored  in  the  capacitor 
is  discharged  over  a  short  duration.  This  parametric  study  indicates 
that  any  improvement  in  the  ionic  conductivity  can  dramatically 
improve  the  supercapacitor  energy  density.  Improving  the  ionic 
conductivity  in  the  electrolyte  is  difficult  and  recent  efforts  are 
to  improve  the  ion  supply  by  electrolyte  additives  [41]  or 
using  transport  [42],  Any  improvements  to  ionic  transport  in  the 
electrolyte  can  drastically  improve  the  energy  stored  in  a 
supercapacitor. 


500 
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Fig.  9.  Current  (during  galvanostatic  charge  and  potentiostatic  discharge) 


for  the  parametric  variation  of  electrode  electronic  conductivity  a. 
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Table  4 

List  of  cases  for  parametric  study  of  k. 


Parameter  Base  case  Test  case-1 


ff(Sm-')  52.1  52.1 

k  (S  nr1)  0.067  0.1 

aC  (F  m-3)  86  x  106  86  x  106 

ts  0.67  0.67 

e,  0.6  0.6 


In  the  previous  two  parametric  sweeps,  we  have  changed  the 
intrinsic  transport  properties  of  the  electrode/electrolyte.  The 
pore  structure  and  connectivity  can  also  influence  the  overall 
ionic  and  electronic  conductivity.  The  effective  transport  property 
is  proportional  to  r\  where  C  is  the  Bruggeman’s  coefficient.  For 
random  packed  spheres,  C  is  1.5  but  it  is  typically  larger  than  1.5 
and  can  be  as  high  as  4.0.  Here  we  vary  the  Bruggeman’s  coef¬ 
ficient  between  1.5— 4.0  to  study  its  effect  on  overall  performance. 
The  higher  values  of  £  implies  higher  tortuosity  and  lower 
effective  transport  coefficients.  We  have  applied  these  Brugge¬ 
man’s  coefficient  for  both  electronic  and  ionic  conductivity  and 
from  the  previous  parametric  study,  the  modification  of  ionic 
conductivity  will  have  the  most  effect  and  Fig.  11  shows  that  the 
impact  to  performance  is  analogous  to  reduction  in  ionic  con¬ 
ductivity  (Fig.  10).  This  suggests  that  any  reduction  in  tortuosity 
of  pore  structure  can  dramatically  improve  the  energy  stored  and 
discharged  and  if  these  structures  have  nearly  extremely  low 
tortuosity  (£  =  1)  as  in  ID  transport,  the  performance  can  be 
further  improved. 

In  the  first  three  parametric  studies,  we  looked  at  the  varia¬ 
tion  of  the  relevant  parameters  controlling  the  transport  in  the 
electrode  and  electrolyte.  The  total  capacity  of  the  super¬ 
capacitors  is  directly  related  to  specific  capacitance,  aC  and  that  is 
the  final  parametric  study  we  have  conducted.  Specific 


capacitance  not  only  depends  on  the  interfacial  surface  area  to 
the  volume  but  also  the  ability  to  store  charge  at  the  double  layer. 
This  parameter  is  highly  dependent  on  the  porosity  of  the  elec¬ 
trode,  the  pore  structure,  the  size  of  the  pores,  the  chemical 
composition  of  the  electrode  and  electrolyte  and  defines  the  total 
capacity  or  energy  that  can  be  stored  in  the  supercapacitor.  In  this 
study,  we  picked  three  values  for  aC  as  shown  in  Table  5.  As 
expected  (Fig.  12),  increasing  the  capacity  increases  the  energy 
stored  in  the  supercapacitor  and  decreasing  the  capacity  as  the 
opposite  effect.  This  is  a  linear  relationship  if  everything  else  is  a 
constant  and  there  are  several  ongoing  efforts  to  increase  the 
specific  capacity  (e.g.,  [14,41]). 

3.3.  2D  and  3D  effects 

3.3.1.  Anisotropy 

Physical  characteristics  such  as  anisotropy  in  various  transport 
properties  (e.g.,  electronic  conductivity)  can  have  significant  effects 
in  terms  of  local  variation  of  the  electrode/electrolyte  potentials, 
current  densities,  species  concentration,  etc.  From  a  system 
perspective,  it  is  important  to  understand  how  that  impacts  the 
overall  cell  performance  and  efficiency.  For  e.g.,  various  meso- 
porous  carbon  materials  show  unusually  high  anisotropic  elec¬ 
tronic  conduction  from  in-plane  conduction  versus  the  through¬ 
thickness  conduction.  We  employ  the  formulation  we  have  devel¬ 
oped  to  these  scenarios  as  we  can  explicitly  account  for  anisotropy 
and  look  at  system  level  performance.  In  Fig.  13,  we  compare  the 
overall  cell  performance  as  a  function  of  time  for  the  isotropic  and 
anisotropic  case.  In  the  anisotropic  case,  we  have  assumed  an 
extreme  case  where  axx  =  52.1  S  nr1  and  ayy  =  azz  -  cr^/lOO.  Even 
in  this  extreme  case,  the  anisotropic  electronic  conductivity  does 
not  have  much  impact  on  the  cell  performance  as  the  transport  is 
predominantly  in  the  ID  x  direction  and  lateral  components  do  not 
affect  the  overall  cell  performance.  In  addition,  as  shown  in  previ¬ 
ous  section,  electronic  conductivity  is  not  the  rate  limiting  property 
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Fig.  10.  Current  (during  galvanostatic  charge  and  potentiostatic  discharge)  versus  time  for  the  parametric  variation  of  electrolyte  ionic  conductivity  ft. 


Fig.  11.  Current  (during  galvanostatic  charge  and  potentiostatic  discharge) ' 


for  the  parametric  variation  of  Bruggeman’s  coefficient 
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Table  5 

List  of  cases  for  parametric  study  of  aC. 

Parameter  Base  case  Test  case-1  Test  case-2  Test  case-3 

(TtSm-1)  52.1  52.1  52.1  52.1 

«(Sm4)  0.067  0.067  0.067  0.067 

aC  (F  m-3)  86  x  106  10  x  106  50  x  106  100  x  106 

es  0.67  0.67  0.67  0.67 

e,  0.6  0.6  0.6  0.6 


and  small  perturbations  in  electronic  conductivity  would  not  have 
much  impact  on  the  overall  cell  performance. 

3.3.2.  Spatial  variability  of  conductivity/void  fraction  in  a  2D 
supercapacitor 

In  this  study  we  introduce  random  variation  around  20%  of  the 
void  fraction  in  the  solid  phase  and  the  corresponding  changes  in 
the  liquid  phase.  We  also  introduce  the  anisotropy  to  understand 
the  interaction  between  spatial  variability  and  anisotropy  and  the 
resulting  gradients  in  the  transverse  direction.  The  Fig.  14  show 
contours  of  variation  of  the  volume  fractions  (with  respect  to  the 
mean  volume  fraction)  of  the  solid  phase. 

Fig.  15  shows  quantitatively  the  change  in  the  solution  profile 
at  the  moment  of  switching  from  charging  to  discharging  process. 
The  error  bars  indicate  the  %  change  in  the  standard  deviation  of 
solution  compared  to  isotropic  case  with  uniform  void  fraction. 
Here  we  note  that  change  is  prominent  during  the  discharging 
process  compared  to  the  charging  process.  This  could  be  because 
of  the  high  gradients  in  the  current  density  introduced  when 
capacitor  is  switched  from  galvanostatic  charging  to  potentio- 
static  discharging  as  can  be  seen  in  Fig.  16.  We  compare  this 
discharging  process  to  the  isotropic  case  and  as  can  be  noted,  the 
2D  Random  and  2D  Random  with  anisotropy  almost  have  similar 
charge  discharge  profiles  indicating  that  even  20%  variability  in 


the  porosity  does  not  have  significant  effect  on  the  overall  cell 
performance  (Fig.  17). 


3.3.3.  Non-planar  supercapacitor 

So  far  we  have  only  shown  the  behavior  of  the  new  formulation 
for  planar  supercapacitor  configurations  where  a  ID  formulation  is 
sufficient  to  model  the  dominant  system  level  processes.  In  this 
section  we  demonstrate  the  inherent  ability  of  the  current  config¬ 
uration  to  model  non-planar  configurations.  For  this  purpose,  we 
consider  an  alternate  configuration  for  the  supercapacitor  where 
the  current  collectors  are  placed  on  the  top  and  bottom  of  the 
negative  and  positive  electrodes  respectively  (shown  in  Fig.  18).  In 
this  configuration  the  ion  transport  is  no  longer  one  dimensional 
and  the  multidimensional  formulation  has  to  be  employed  to 
simulate  these  non-planar  configurations. 

In  Figs.  19  and  20  we  compare  the  spatial  distribution  of  cur¬ 
rent  immediately  after  we  switch  to  potentiostatic  discharge 
between  the  planar  configuration  and  non-planar  configuration. 
In  the  planar  configuration,  the  current  density  contour  surfaces 
are  not  perfectly  planar  close  to  the  current  collector  but  quickly 
become  planar  near  the  separator.  This  indicates  that  the  pre¬ 
dominant  ion  and  electron  transport  becomes  essentially  ID  even 
at  these  high  discharge  rates  after  the  switch.  In  comparison  the 
non-planar  case  has  2d  curved  surfaces  over  the  entire  electrode 
region  with  some  straightening  over  the  narrow  separator.  As  one 
can  see  higher  current  densities  closer  to  separator  as  compared 
to  far  off  regions.  Visually  one  can  also  see  better  utilization  of 
the  electrode  in  the  2D  configuration  (similar  to  multidimen¬ 
sional  architectures  in  batteries).  As  in  batteries,  novel  2D  and  3D 
configurations  offer  ways  to  reduce  the  transport  distance  be¬ 
tween  the  positive  and  negative  electrodes  while  ensuring  better 
volumetric  utilization  of  the  electrodes.  In  future,  we  will 
perform  optimization  studies  to  arrive  at  novel  configurations 
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Fig.  12.  Current  (during  galvanostatic  charge  and  potentiostatic  discharge)  versus  time  for  the  parametric  variation  of  specific  capacitance  (aC). 


Fig.  13.  Current  (during  galvanostatic  charge  and  potentiostatic  discharge) ' 


:  for  the  anisotropic  vs.  isotropic  electronic  conductivity  of  the  electrode. 
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Fig.  14.  Random  Variations  in  solid  phase  void  fractions  normalize  with  respect  to  the  mean  volume  fraction. 


and  it  is  always  a  challenge  how  one  can  fabricate  these  complex 
supercapacitor  assemblies. 

4.  Conclusions 

A  generalized  three-dimensional  formulation  for  the  electric 
double  layer  supercapacitor  has  been  proposed,  implemented,  and 


tested  in  this  paper.  Simulation  results  are  in  very  good  agreement 
with  the  published  experimental  data.  The  three-dimensional 
nature  of  the  charge  transport  and  current  densities  at  small 
scales  due  to  spatial  variations  in  porous  electrodes  has  been 
revealed.  These  spatial  variations  do  not  have  significant  variations 
in  the  global  response  of  the  cell  as  the  predominant  diffusion 
processes  smear  the  local  effects.  For  standard  supercapacitor  with 


Fig.  15.  Variation  i 
and  anisotropy. 


i  spatial  variability 


:  al.  /  Journal  of  Power  Sources  256  (2014)  369—382 


Fig.  16.  Variation  in  current  density  at  various  X  locations  during  the  initial  phase  of  the  discharging  process  along  with  %  deviation  showing  in-plane  variation  of  the  current 
density  due  to  spatial  variability  and  anisotropy. 
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Fig.  17.  Comparison  of  random  variation  of  solid  and  liquid  state  void  fractions. 
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Fig.  18.  Schematic  of  the  non-planar  supercapacitor. 


planar  electrodes,  the  results  of  the  3D  model  are  analogous  to  ID 
models  but  the  generality  of  the  3D  modeling  framework  can 
capture  the  local  spatiotemporal  variations  and  serve  as  a  direct 
link  to  the  material  properties  at  the  mesoscale.  In  addition,  these 
models  also  provide  the  capability  to  design  new  cell  configura¬ 
tions  that  provide  additional  flexibility  to  minimize  the  transport 
distance  between  the  electrodes,  packing  more  electrode  material 
and  thereby  improving  the  power/energy  density  of  the  super¬ 
capacitors.  In  addition,  the  degradation  is  likely  to  be  a  local 
phenomena  and  having  a  3D  local  model  makes  it  amenable  to 
incorporate  mechanistic  degradation  mechanisms.  It  is  also  found 
from  the  sensitivity  analysis  that  the  supercapacitor  performance 
is  most  sensitive  to  specific  capacitance,  liquid-phase  ionic 
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Fig.  19.  Spatial  distribution  of  the  current  density  for  parallel  plate  supercapacitor  configuration. 


Fig.  20.  Spatial  distribution  of  the  current  density  for  non-planar  supercapacitor  configuration. 


conductivity,  and  Bruggeman’s  coefficient.  This  indicates  that 
system  level  improvements  can  be  made  by  designing  materials  to 
have  high  specific  capacitance,  improved  ionic  conductivity  and 
very  low  tortuous  electrodes. 
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